A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11 (09.01.1; 09.13.2; 11.01.1; 11.13.1; 12.03.3; 13.19.3) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



ON 
ON 



Molecular abundances in the Magellanic Clouds 
II. Deuterated species in the LMC 

Y.-N. Chini 2, C. Henkel^, T.J. Millar'', J.B. Whiteoak^ and R. Mauersberger^^ 

^ Institute of Astronomy and Astrophysics, Academia Sinica, P.O.Box 1-87, Nankang, 115 Taipei, Taiwan 
Radioastronomisclies Institut der Universitat Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
Max-PIanck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

* Department of Physics, UMIST, P O Box 88, Manchester M60 IQD, United Kingdom 

^ Paul Wild Observatory, Australia Telescope National Facility, CSIRO, Locked Bag 194, Narrabri NSW 2390, Australia 
^ Steward Observatory, The University of Arizona, Tucson, AZ 85721, U.S.A. 



=3 

1—5 

o 



> 
cn 
in 
o 
i> 
o 
\o 

I 

o 



Received date ; accepted date 



Abstract. The first definite discoveries of extragalactic 
deuterium are reported. DCO+ has been detected in three 
and DCN has been measured in one star-forming re- 
gion of the Large Magellanic Cloud (LMC). While the 
HCO+ /DCO+ abundance ratios are found to be 19 ± 3, 24 
± 4, and 67 ± 18 for N113, N44BC and N159HW, respec- 
tively, a HCN/DCN abundance ratio of 23 ± 5 is obtained 
for N113. These results are consistent with a gas temper- 
ature of about 20 K and a D/H ratio of about 1.5 x 10^^, 
consistent with that observed in the Galaxy. If the cloud 
temperature is closer to 30 K, then a D/H ratio is required 
to be up to an order of magnitude larger. Because this ra- 
tio provides a lower limit to the primordial D/H ratio, it 
indicates that the baryon mass density alone is unable to 
close the universe. 

Key words: ISM: abundances - ISM: molecules ~ Galax- 
ies: abundances - Magellanic Clouds - Cosmology: obser- 
vations - Radio lines: ISM 



1. Introduction 

The primordial ratio of the number of deuterium to hy- 
drogen nuclei (D/H) created in big bang nucleosynthesis is 
the most sensitive measure of the cosmological baryon-to- 
photon ratio, 77, and the cosmological density of baryons, 
Qb- In the interstellar medium of our Galaxy D/H = 
1.6 ±0.3 X 10~^ (Linsky et al. 1995), which places a strict 
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* Based on observations with the Swedish-ESO Submillime- 
ter Telescope (SEST) at the European Southern Observatory 
(ESO), La Silla, Chile 



lower limit on the primordial abundance, because stars re- 
duce the proportion of D in the ISM. On the other hand, 
quasar absorption line systems should provide suitable 
targets to determine primordial D abundances because 
they can sample metal-poor gas at early epochs where 
the destruction of D should be negligible. Carswell et al. 
(1994) and Songaila et al. (1994) obtained a high D/H 
w 2.5 X 10"'* ratio from a possible deuterium absorption 
feature toward the quasar Q0014+813. More recent ob- 
servations by Rugers & Hogan (1996) and Tytler et al. 
(1996) do not provide clear evidence in favor or against 
high D/H ratios at large redshifts and it is still difficult to 
rule out entirely the possibility that the D feature is con- 
taminated by H having a velocity different from that of 
the principal hydrogen lines. Therefore, a direct measure- 
ment of atomic deuterium toward high redshift quasars 
is a potentially powerful, though not yet fully explored 
means of studying D/H. 

Irrespective of the validity of a high or low cosmo- 
logical D/H value, the ratio appears to be small enough 
to preclude any easy detection of deuterated molecu- 
lar species in interstellar space. Nevertheless, the zero- 
point energy difference between hydrogen- and deuterium- 
bearing molecules means that deuterated species can be- 
come heavily fractionated in cool environments. Models 
of fractionation (Brown & Rice 1986; Millar et al. 1989) 
indicate that the underlying D/H ratio is consistent with 
that suggested by Linsky et al. (1995), implying that the 
baryon mass density is not large enough to close the uni- 
verse. Hence measurements of deuterated species have be- 
come an important tool to study the physical properties 
and the chemical and dynamical history of Galactic molec- 
ular clouds (e.g. Wootten 1987; MiUar et al. 1989). In a 
previous study Mauersberger et al. (1995) presented sensi- 
tive upper limits to the DCN/HCN ratio in NGC 253 and 



Y.-N. Chin, C. Henkel, T.J. Millar, J.B. Whiteoak et al.: Molecular abundances in the Magellanic Clouds 



MmK] N113 

10 
5 


-5 

200 
100 




N113 




200 220 240 260 
Velocity (km/s) 



200 220 240 260 
Velocity (km/s) 



,^^[^ K] N44BC 

hi I I I I I I I 

-DCO (2-'l) 



^MBLmK] 



N159HW 



20 
10 


400 

200 


15 
10 

5 


-5 




260 280 300 320 
Velocity (km/s) 



220 240 260 
Velocity (km/s) 



Fig. 1. The observed spectra of selected hydrogen-, deuterium-, and ^^C-bearing molecules toward N113 (aigso = 5'^13™38'!7, 
51950 = -69°25'57"), N44BC (aigso = 5''22"^10=6, ^igso = -68°00'32"), and N159HW (aigso = 5''40'"04':4, Sigso = -69°46'54") 
in the LMC 



IC342, which are consistent with the warm, metal rich 
environment of this starburst nucleus. 



2. Observations 

The data were taken in January and March 1996, using the 
15-m Swedish-ESO Submillimetre Telescope (SEST) at La 
Silla, Chile. Two SIS receivers, one at A = 3 and one at 
2 mm, yielded overall system temperatures, including sky 
noise, of order Tgys — 250 K on a main beam brightness 
temperature (Tmb) scale. The backend was an acousto- 
optical spectrometer which was split into 2 x 1000 con- 
tiguous channels for simultaneous 3 and 2 mm observa- 
tions. The channel separation of 43 kHz corresponds to 
0.09 and 0.15 km s"^ at 145 and 88 GHz, respectively. The 
corresponding antenna beamwidths were 35" and 55" at 
the observed line frequencies taken from Lovas (1992). 

The observations were carried out in a dual beam- 
switching mode (switching frequency 6 Hz) with a beam 
throw of 11'40" in azimuth. The on-source integration 
time of each spectrum varied from 8 minutes for HCO+ 
to 18 hours for H^^CN toward N113. All spectral intensi- 
ties obtained were converted to a Tmb scale, correcting for 
main beam efficiencies of 0.74 at 85-100 and 0.67 at 130- 
150 GHz (Dr. L.B.G. Knee, priv. comm.). The pointing 
accuracy, obtained from measurements of the SiO maser 
sources RDor and UMen, was better than 10". 



3. Results 

Toward the prominent star- forming region N113 we have 
detected HCO"*", HCN, and their rare isotopic species 
H13C0+, DCO+, Hi^CN, and DCN. HCO+, Il^^CO+, 
and DCO+are also observed toward N159HW (position 
defined by Hunt & Whiteoak 1994) and N44BC. The spec- 
tra and line parameters obtained from Gaussian fits are 
displayed in Fig. and Table |l|. 

Assuming optically thin line emission, we can derive 
from the line intensity ratios of the ^^C- and deuterium- 
bearing species the ratio of column densities. We find 
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{N: column density; B: rotational constant; gi and g„: 
statistical weights of lower and upper states; T: excitation 
temperature; I: integrated line intensity; indices 1 and 2 
refer to the deuterium-bearing and to the ^^C-bearing or 
main isotopic species, respectively). Both HCO"*" and HCN 
have large electric dipole moments (of order 3 Debye) and 
are excited well beyond 2.7 K only at high densities. At 
n(H2) = 10^ cm^"^ radiative transfer calculations with 
a Large Velocity Gradient (LVG) model, spherical cloud 
geometry, and collision cross sections from Green & Thad- 
deus (1974) yield excitation temperatures of order 4K 
for Tkin = 20-40 K in the optically thin limit, both for 
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Table 1. Parameters of the Observed Molecular Lines 
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a) The hyperfine components of DCN and H^^CN line cannot be resolved. 

b) The three HCN hyperfine transitions {F=l-1, F=2-l, _F=0-1) have been resolved by a Gaussian fit. While Tmb values for 
each component are given, the total integrated line intensity refers to the entire line. 

c) During the Gaussian fit value of Avi/2 has been fixed by using that determined from the H^^CO^(l-O) transition. 

a less 'chemically' processed environment than that pro- 
vided by the Galaxy. Gensheimer et al. (1996) observed 
a high deuteration of water toward a number of Galac- 
tic hot cores despite a high kinetic temperature. While 
this can be explained in terms of a recent evaporation of 
grain mantles, this effect probably plays no role toward 
our LMC sample since we are looking on a much larger 
region which may much better reflect the average condi- 
tions of the interstellar gas than the highly biassed sam- 
ple by Gensheimer et al. (1996). Due to the lower metal- 
licity in the LMC, we expect grain surface reactions to 
be of lesser importance than for the Galactic interstel- 
lar medium. The process of chemical fractionation has, 
however, to be accounted for quantitatively. To this end, 
we have re-visited earlier models of deuterium chemistry 
(Millar et al. 1989; Howe & Millar 1992) with conditions 
appropriate for clouds in the LMC (Paper I) and with an 
elemental D to H ratio of 1.5 x 10~^. We have updated 
the chemistry to include the fast dissociative recombina- 
tion of (and H2D+) and more accurate rate coeffi- 
cients for neutral-neutral reactions. Although the general 
trends of fractionation can be understood, especially for 
simple species like HCO+ detailed models are necessary 
for quantitative comparison with the observations (Millar 
et al. 1989). 



the J— 1-0 transition of H^^CN and for the 2-1 transi- 
tion of DCN. These results also hold for HCO+. While 
smaller densities would lead to excitation temperatures 
almost undistinguishable from those of the microwave 
background (and thus to negligible emission), gas at even 
higher densities is likely not common enough to contribute 
significantly to the observed emission. We thus find, ne- 
glecting differences in beam sizes, for HCO+ 

^^(DCO+) ^ /(DCO+) 

Ar(Hi3C0+) ^ /(H13C0+) ^ ' 
and for HCN 

A^(DCN) ^ /(DCN) 

A^(Hi3CN) ^ /(H13CN)' ^ ' 

With 12C/13C « 50 (Johansson et al. 1994; this is con- 
sistent with the line intensity ratios dereived from Ta- 
ble|) we then obtain iV(HCO+)/iV(DCO+) « 19 ± 
3 and 7V(HCN)/7V(DCN) « 23 ± 5 for N113, while 
7V(HC0+)/iV(DC0+) « 24 ± 4 and 67 ± 18 for N44BC 
and N159HW, respectively. 

4. Discussion 

The detection of DCO+ and DCN in the LMC is a use- 
ful tool to study D/H on a larger spatial scale and in 
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Table 2. Calculated abundance ratios 



Temperature 
Time (yr) 


10 K 




20 K 






30 K 




i?(H+) i?(HCO+) i?(HCN) 


RiHt) 


i?(HCO+) 


i?(HCN) 


R{U+) 


i?(HCO+) 


i?(HCN) 


10^ 


3.8 12.8 24.9 


8.3 


23.7 


61.1 


167 


309 


246 


3.2 X 10^ 


3.7 10.3 14.9 


8.2 


17.7 


32.5 


166 


257 


253 


10^ 


4.3 9.3 20.0 


8.6 


14.8 


37.5 


166 


251 


476 



a) The ratio i?(HX) is defined as the abundance ratio of HX to DX, and is tabulated for HX = H^, HCO"*" and HCN 



Because deuterium fractionation depends on small dif- 
ferences in zero-point energies, the process is sensitive to 
the kinetic temperature, Tkim of the molecular gas. Multi- 
transition studies of CO in N113 have not been made so 
that Tkin is unknown, but a value of around 20 K is ex- 
pected since the cloud is associated with an H ii region 
(for detailed modelling, see Lequeux et al. 1994). To cover 
the range of possibilities, we present the results of time- 
dependent chemical kinetic calculations at 10, 20 and 30 K 
for three model times in Table||. 

At low temperatures there are two main routes to D 
fractionation: 



H 



HD 

+ D 



H2D+ + H2 
H2D+ + H 



(4) 
(5) 



where we have calculated the rate coefficient, fc_4, for the 
back reaction of H2D+ with H2 using the data presented 
by Sidhu et al. (1992). This rate coefficient is very sensitive 
to temperature and has values of 2.1 x 10~^^, 2.4 x 10"^'^ 
and 8.3 x lO'^^ cm^s-^ at 10, 20 and 30 K, respectively. 
For r <; 25 K, the back reaction is fast and i?(H;^) is large, 
« 165 at 30 K. DCO+ is formed through reactions involv- 
ing H2D"'"and D atoms, with these species contributing 
64% (36%) at 10 K, 54% (46%) at 20 K and 44% (56%) at 
30 K, to the formation rate of DCO"*". 

Our results show that the deuterium fractionation is 
larger than that observed in N113 for HCO"'" and HCN 
for r < 20K and smaller for T = 30 K; the observed re- 
sult that fractionation is larger in DCO"^ than DCN is 
reproduced. 

Our results are consistent with T ~ 20K for a D/H 
ratio of 1.5 x lO^'^. However, the temperatures of the 
LMC clouds are not, as yet, well constrained. The ob- 
served abundance ratios could be reproduced in warmer 
clouds if the D/H ratio is larger than our adopted value. 
At a cloud temperature of 30 K, one would need a D/H 
ratio about 10 times larger, i. e. on the order of 1.5 x 10~*. 
If there temperature is 20-30 K as suggested by Lequeux 
et al. (1994) for the cloud cores (for HCN, see also Chin et 
al. 1996), then it appears that the D/H ratio in the LMC is 
similar to or larger than that in our Galaxy. This would be 
consistent with the LMC gas being less processed through 
stars than that in our Galaxy. Thus, the evidence from 



the LMC is consistent with an open universe as long as 
the cosmological constant is small. 

Acknowledgements. We like to thank Drs. J.L. Linsky and T.L. 
Wilson for useful comments. TJM is supported by a grant 
from PPARC. RM was supported by a Heisenberg fellowship 
of Deutsche Forschungsgemeinschaft. 

References 

Brown, R.D., Rice, E.H.N. , 1986, MNRAS 223, 429 
Carswell, R.F., Ranch, M., Weymann, R.J., Cooke, A.J., Well, 

J.K., 1994, MNRAS 268, LI 
Chin, Y.- N., Henkel, C, W; hiteoak, J.B., et al, 1996, A&A, in 



press. 



astro-ph/9606081 (Paper I) 



Gensheimer, P., Mauersberger, R., Wilson, T.L. 1996, A&A, 
in press 

Green, S., Thaddeus, P., 1974, ApJ 191, 653 
Howe, D.A., Millar, T.J., 1992, MNRAS 262, 868 
Hunt, M.R., Whiteoak, J.B., 1994, Proc. Astron. Soc. Aust. 
11, 68 

Johansson, L.E.B., Olofsson, H., Hjalmarson, A., Gredel, R., 

Black, J.H., 1994, A&A 291, 89 
Lequeux, J., Le Bourlot, J., Pineau de Forets, G., et al., 1994, 

A&A 292, 371 

Linsky, J.L., Diplas, A., Wood, B.E., et al., 1995, ApJ 451, 335 
Lovas, F.J., 1992, J. Phys. Chem. Ref. Data, 21, 181 
Mauersberger, R., Henkel, C, Chin, Y.-N., 1995, A&A 294, 23 
Millar, T.J., Bennett, A., Herbst, E., 1989, ApJ 340, 906 
Rugers, M., Hogan, C.J., 1996, ApJ 459, LI 
Sidhu, K.S., Miller, S., Tennyson, J., 1992, A&A 255, 453 
Songaila, A., Cowie, L.L., Rugers, M., 1994, Nat 368, 599 
Tytler, D., Fan, X.-M., Buries, S., 1996, Nat 381, 207 
Watson, W.D., 1976, Rev. Mod. Phys. 48, 513 
Wootten, A., 1987, lAU Symp. 120, Astrochemistry, eds. 
M.S. Vardya and S.P. Tarafdar, D. Reidel Publishing Com- 
pany, Dordrecht, p. 311 



This article was processed by the author using Springer- Verlag 
X^TgX A&A style file L-AA version 3. 



